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ABSTRACT 


Since the beginning of its routine science operation in March 
2015, the NASA SMAP observatory has been returning 
interference-mitigated brightness temperature observations 
at L-band (1.41 GHz) frequency from space. The resulting 
data enable frequent global mapping of soil moisture with a 
retrieval uncertainty below 0.040 m?/m? at a 36 km spatial 
scale. 

This paper describes the development and validation of 
an enhanced version of the current standard soil moisture 
product. Compared with the standard product that is posted 
on a 36 km grid, the new enhanced product is posted on a 9 
km grid. Derived from the same time-ordered brightness 
temperature observations that feed the current standard 
passive soil moisture product, the enhanced passive soil 
moisture product leverages on the Backus-Gilbert optimal 
interpolation technique that more fully utilizes the additional 
information from the original radiometer observations to 
achieve global mapping of soil moisture with enhanced 
clarity. The resulting enhanced soil moisture product was 
assessed using long-term in situ soil moisture observations 
from core validation sites located in diverse biomes and was 
found to exhibit an average retrieval uncertainty below 
0.040 m?/m?. As of December 2016, the enhanced soil 
moisture product has been made available to the public from 


the NASA Distributed Active Archive Center at the National 
Snow and Ice Data Center. 


1. INTRODUCTION 


Since the beginning of its routine science operation in March 
2015, the NASA Soil Moisture Active Passive (SMAP) 
observatory [1] has been returning interference-mitigated 
brightness temperature (Tp) observations at L-band 
frequency (1.41 GHz) from space. The resulting data enable 
the production of the standard SMAP Passive Soil Moisture 
Product (L2_SM_P) [2] that provides frequent global 
mapping of soil moisture with a retrieval uncertainty below 
0.040 m?/m? at a 36 km spatial scale [3]. 

Given the scanning geometry of the radiometer onboard 
the SMAP observatory, the time-ordered Tz observations in 
the standard Level 1 Brightness Temperature Product 
(L1B_TB) [4] are oversampled (i.e., exceeding the Nyquist 
sampling rate) in the along-scan direction. By applying the 
Backus-Gilbert (BG) optimal interpolation technique [5][6] 
to LIB_TB, additional radiometric information that was not 
readily available in the current standard Level 1 Gridded 
Brightness Temperature Product (L1C_TB) [7] can be more 
fully utilized for the development of higher-level 
geophysical products, such as the SMAP Enhanced Passive 
Soil Moisture Product (L2_SM_P_E). 


The development of L2_SM_P_E begins with the 
application of the BG optimal interpolation technique to the 
time-ordered Tg observations in LIB_TB. This technique 
makes use of the antenna gain pattern information specific to 
the 6-meter reflector assembly onboard the observatory to 
achieve optimal data interpolation on a fine grid from the 
original data. The interpolation is considered optimal in the 
sense that it aims at producing interpolated values closest to 
what would have been measured had the instrument actually 
made its measurements at the interpolation points [5]. A 
comparison between the standard L1C_TB on a 36 km grid 
and the SMAP Level 1 Enhanced Gridded Brightness 
Temperature Product (L1IC_TB_E) on a 9 km is shown in 
Fig. 1. 


(Eee 
220 224 228 232 236 240 244 248 252 256 260 264 268 272 276 280 284 288 292 296 300 


(b) 


Figure 1: SMAP horizontally polarized TB observations 
during December 15-17, 2016 over the Amazon basin. The 
application of BG interpolation in L1C_TB_E (Fig. 1a) 
results in enhanced spatial features of the scene that are often 
concealed or not immediately obvious in the standard 
L1C_TB product (Fig. 1b). 


It is clear from the comparison that the enhanced 
L1C_TB_E product (Fig. la) reveals spatial features of the 
scene that are concealed or not immediately obvious in the 
current standard L1C_TB product (Fig. 1b). Overall, the 
L1C_TB_E image presents a less pixelated representation of 
the original data, offering a higher acuity and contrast 
compared with the LIC_TB image. Because the 
information in the along-scan oversampled measurements in 
the original time-ordered Tg data is more completely 
captured in LIC_TB_E, the product provides a more 
optimal representation of the original data for the production 
of L2_SM_P_E. 


2. PRODUCT DEVELOPMENT 


The development of the enhanced L2_SM_P_E follows that 
of the standard L2_SM_P [3]. A _ flowchart of the 
L2_SM_P_E processing is shown in Fig. 2. 
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Figure 2: Processing flow of the L2_SM_P_E processor. The 
processor uses the L1C_TB_E product as the primary input, 
along with other ancillary data on finer grid resolutions, to 
retrieve soil moisture. 


Both products share the same baseline soil moisture 
retrieval algorithm — the single channel algorithm that uses 
the vertically polarized Tg channel (SCA-V) and the same 
set of ancillary data. The fore- and aft-look Tg observations 
in L1C_TB_E are first combined to provide the primary 
input data to the L2_SM_P_E processor. External static and 
dynamic ancillary data preprocessed on finer grid 
resolutions are then brought into the processing to evaluate 
the feasibility and subsequent estimated quality of the 
retrieval. | When surface conditions favorable to soil 
moisture retrieval are identified at a given grid cell, retrieval 
is performed. Corrections for water contamination, surface 
roughness, effective soil temperature, and vegetation water 
content are applied to produce the final output soil moisture 
retrieval field on the same 9 km EASE Grid 2.0 global 
projection [8][9] as the input LIC_TB_E. The resulting 
data contents of the L2_SM_P_E output files are described 
in the Product Specification Document [10] of the product. 

Compared with the standard L2_SM_P that is posted on 
a 36 km grid, the enhanced L2_SM_P_E is posted on a 9 km 
grid. As such, the processing for L2_SM_P_E is (36/9)? or 
16 times more intensive than that for L2_SM_P. Despite 
this increase in computation time, L2_SM_P_E processing 
remains relatively swift on the SMAP Science Data System 
(SDS). 


3. PRODUCT VALIDATION 


The retrieval performance of L2_SM_P_E was assessed 
using the same set of validation methodologies described in 
[3]: Long-term ground truth consisting of properly scaled 
aggregations of in situ soil moisture observations from core 
validation sites serves as a primary way for product 
assessment. Agreement between the L2_SM_P_E soil 
moisture estimates and in situ data over space and time are 
reported in time-series correlation coefficient, bias, root- 
mean-square difference (RMSD), and unbiased root-mean- 
square difference (ubRMSD). Together, these metrics 
provide a more comprehensive description of product 
performance than any one alone. Among these metrics, the 
ubRMSD (in units of m?/m*) is the metric adopted by SMAP 
for stating product accuracy across all Level 2 through Level 
4 soil moisture products. The SMAP Level | mission 
requirement for the active/passive soil moisture product 
accuracy is ubRMSD = 0.040 m*/m3 or better. The same 
accuracy target was also adopted for the standard L2_SM_P. 

In situ soil moisture observations from a total of 15 core 
validation sites were used to assess the retrieval performance 
of L2_SM_P_E at a spatial scale similar to the effective 
field-of-view (EFOV) of the SMAP radiometer. 
Performance metrics were evaluated using 19 months (April 
1, 2015 through October 31, 2016) of in situ soil moisture 
observations from core validation sites and _ the 
corresponding L2_SM_P_E granules from both 6:00 am 
descending passes and 6:00 pm ascending passes. Tables 1 
and 2 provide a list of the core validation sites used in this 
assessment and the corresponding performance metrics, 
tabulated separately for the 6:00 am / 6:00 pm L2_SM_P_E 
data. 


Table 1: Assessment of 6:00 am descending L2_SM_P_E 
based on a 19-month matchup with in situ soil moisture 
observations from 15 core validation sites. The units for 
ubRMSD, RMSD, and bias are m?/m’. 


Sites ubRMSD Bias RMSD Corr 
Reynolds Creek 0.040 -0.023 0.046 0.598 
Walnut Gulch 0.024 0.011 0.026 0.813 
TxSON 0.032 -0.015 0.036 0.921 
Fort Cobb 0.028 -0.056 0.062 0.883 
Little Washita 0.022 -0.027 0.035 0.920 
South Fork 0.054 -0.062 0.082 0.646 
Little River 0.028 0.087 0.092 0.887 
Kenaston 0.022 -0.040 0.046 0.854 
Carman 0.056 -0.088 0.104 0.611 
Monte Buey 0.051 -0.020 0.055 0.840 
REMEDHUS 0.042 -0.007 0.042 0.872 
Twente 0.056 0.013 0.057 0.885 
HOBE 0.036 -0.009 0.037 0.863 
MAHASRI 0.036 -0.006 0.037 0.728 
Yanco 0.043 0.020 0.048 0.964 
L2_SM_ PE 0.038 -0.015 0.054 0.819 
L2_SM_P 0.037 -0.014 0.052 0.822 


A few observations can be drawn from the performance 
metrics listed in Tables 1 and 2. First, the enhanced 


L2_SM_P_E and the standard L2_SM_P showed great 
similarity in all performance metrics; both demonstrated an 
average ubRMSD below 0.040 m*/m? and thus met the 
SMAP Level 1 mission requirement for soil moisture 
retrieval accuracy. The performance of the 6:00 am data is 
slightly better than that of the 6:00 pm data. This is within 
expectations because the current correction in geophysical 
inversion for effective soil temperature described in [3] had 
been optimized primarily for the 6:00 am data only. Despite 
this limitation that led to slight degradation in retrieval 
accuracy the performance metrics for both 6:00 am data and 
6:00 pm data are nonetheless similar. 


Table 2: Assessment of 6:00 pm ascending L2_SM_P_E 
based on a 19-month matchup with in situ soil moisture 
observations from 15 core validation sites. The units for 
ubRMSD, RMSD, and bias are m?/m’. 


Sites ubRMSD Bias RMSD Corr 
Reynolds Creek 0.042 -0.042 0.059 0.651 
Walnut Gulch 0.029 -0.019 0.034 0.676 
TxSON 0.028 -0.018 0.034 0.929 
Fort Cobb 0.035 -0.069 0.077 0.846 
Little Washita 0.026 -0.032 0.041 0.910 
South Fork 0.045 -0.087 0.098 0.764 
Little River 0.029 0.078 0.083 0.872 
Kenaston 0.027 -0.051 0.057 0.828 
Carman 0.049 -0.109 0.120 0.594 
Monte Buey 0.052 -0.019 0.056 0.874 
REMEDHUS 0.045 -0.018 0.048 0.857 
Twente 0.052 0.001 0.052 0.903 
HOBE 0.042 -0.013 0.044 0.844 
MAHASRI 0.038 -0.018 0.042 0.700 
Yanco 0.053 0.011 0.054 0.966 
L2_SM_P_E 0.039 -0.027 0.060 0.814 
L2_SM_P 0.039 -0.028 0.061 0.795 


Overall, the development of L2_SM_P_E results in a 9 
km global soil moisture product that benefits from the 
enhanced information available in the oversampled time- 
ordered Tg observations in the original data, while retaining 
decent retrieval performance (ubRMSD below 0.040 m?/m?) 
as the standard L2_SM_P. 


4. CONCLUSION 


In this paper, the development and validation of the SMAP 
Enhanced Passive Soil Moisture Product were discussed. 
Posted on a 9 km grid on EASE Grid 2.0 global projection, 
the new product leverages on the Backus-Gilbert optimal 
interpolation technique that more fully utilizes the additional 
information from the original radiometer observations to 
achieve global mapping of soil moisture with enhanced 
clarity. Assessment of the enhanced product based on 19 
months of in situ soil moisture observations from core 
validation sites showed that the enhanced product exhibits a 
level of retrieval performance (ubRMSD below 0.040 
m?/m?) that met the SMAP Level 1 mission requirement for 
soil moisture retrieval accuracy. As of December 2016, the 


enhanced soil moisture product has been made available to 
the public from the NASA Distributed Active Archive 
Center (DAAC) at the National Snow and Ice Data Center 
(NSIDC). 
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